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Abstract. 

Exact loop-variables formulation of pure gauge lattice QCD3 is derived from the Wilson 
version of the model. The observation is made that the resulting model is two-dimensional. 
This significant feature is shown to be a unique property of the gauge field. The model is 
defined on the infinite genus surface which covers regularly the original three-dimensional 
lattice. Similar transformation applied to the principal chiral field model in two and three 
dimensions for comparison with QCD. 
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QCD is believed to be a genuine model of the strong interaction. This belief is mainly 
based on the discovery of asymptotic freedom Q and on the subsequent success of the pertur- 
bative description of ultra-violet physics. However, the main mystery of the strong interac- 
tions, its infra-red behavior and, in particular, the confinement of quarks, remains unexplained 
because of the lack of a non-perturbative solution. A number of 'mechanisms' of confinement 
have been suggested. However, no direct link has been unambiguously established between 
the model itself and any particular 'mechanism'. Nevertheless, the model is still analyzable 
by direct methods. 

The unique, in many respects, non-perturbative formulation of QCD is the Wilson lattice 
model The partition function of the pure gauge model in D dimensions has the form 

z= [ llwlldu^frv^sfaYlUt) , (1) 

i v ie P 

where A D is the bare coupling constant, I and p denote links and plaquettes of D-dimensional 
lattice, Ui is the unitary matrix (U(N) or SU(N)) attached to l-th link. In all the products 
the link variables are ordered. 

The gauge-invariant ^-function is 

s{u P ,Y[Ut) =£x.(£4)xr(n^) . ( 2 ) 

I r I 

where r is an irreducible representation with character Xr(U). 

The only non-zero observables in this model are invariant ordered products of link vari- 
ables along closed loops, such as 

W (C) = <-U J] U t ) (3) 

zee 

(C is the closed contour). It has been understood that the loop variables are relevant to 
solution of the quark confinement problem [||, ||. In particular, as argued in Q, the simplest 
criterion of confinement is the area-law behavior of (||): W(C) ~ exp(— a A), where A is the 
area of the 'minimal' surface bounded by C and a is the positive parameter (string tension) 
to be computed. On a lattice level, A is the number of plaquettes of a disk bounded by 
C. Since the loop variables are only non-zero observables, it is clear that the model can be 
formulated in terms of these variables only. There is a number of works devoted to the loop 
dynamics (see review Q and references therein). However, this has not led to a realization 
of confinement. In contrast, one can argue that if the model is a free theory in terms of 
the plaquette variables, then the area law is straightforward. Such realization requires direct 
reformulation of the partition function (|l|) to the loop (plaquette) variables. The latter is the 
main purpose of present paperQ. 

To rewrite partition function (|l]) in terms of the plaquette variables only, we have to 
perform integration over all link variables. 

First, we replace @ by 

S(l, U} U U l)=J2 d rXr(ul J] , (4) 
l_ r I 

x See note added to the end of this paper. 
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Figure 1: Choice of A- and B- vertices and the hexagonal section of a cube. 



where d r = Xr{I) is the dimension of representation. Such replacement does not change (||) 
due to invariance of the action with respect to similarity transformations U p — > QU P & and 
due to the formula 

J cM X r(nArtB) = l Xr (A)xr(B) . (5) 

However, there is a considerable difference in applications. In the case of definition (||), the 
O-integrals are completely decoupled, while in the case of (^), they play an important role. 
The reason we choose the definition (||), though at a first glance it seems an unnecessary 
complication, will soon become clear. Actually, we will use the freedom to choose the point 
of insertion of III into the product of t/j's in (Q), which means that we mark the point where 
the product around the plaquette starts (and ends). 

Now we change the order of integration in (jl]) and consider the (topological) integral over 
link variables 

/ n ^ip(^n^) > w 

/ p ie P 

where A4 is the number of /c-simplexes: Ao is the number of vertices, A/i is the number of 
links etc. The number of independent integrals is Ai — A/"o due to the gauge invariance. 

In this paper, we only consider the 3d caseQ. Due to the Euler theorem for the regular 
infinite 3d lattice, 

AA -M+AA 2 -AA 3 = 0, (7) 

the integral @ must be proportional to the product of A/3 ^-functions of UpS. 

Thus we have to find the set of A/3 constraints (one for each cube) self-consistent with the 
original definition of the model. Let us choose two vertices, A and B, of any cube, that are 
maximally remote from each other (Fig.|l|), and imply the choice is valid for the neighboring 
cubes, i.e. the ^-vertex is the same for all four cubes it belongs to, and similarly for B. Then, 
for the rest of the lattice, all A- and S-vertices are uniquely defined. The boundary of any 
plaquette passes through either an A- or a S-vertex and we call the corresponding plaquette 
either a pa- or a p^-plaquette. Now we use our freedom in placing the f/ p -matrices into the 
product around the plaquette^ so that the product of links around the pA-plaquettes starts 

2 The 2d model is exactly solvable || and its general solution has been found in The 4d case will be 
considered in the next paper Q. 
3 See comments to (W). 
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Figure 2: (a) The 8-cube's fragment of the lattice, ^-vertices are the same as in Fig.l. Hexagons (sections 
of cubes) form the 2d lattice (not all the lines are shown in order to not disturb the view), (b) The fragment 
of the resulting (smoothed) surface (corresponds to 32 cubes of original lattice) . 



at the vertex A, and similarly for B. Then, the invariant ordered product of three U PA s of 
each cube is equal to one of t7 ps 's of the same cube, i.e., 



Xr(l[U pA )=Xr(l[U 1 

PA PB 



pb ; ' 



and both are equal to the invariant ordered product of t/j's, where Us form the closed loop 
separating A and B on the cube surface (double line in Fig.Q). Then, the sought-for constraint 
is expressed by the formula: 



PA Pb r p A p B 



(8) 



where the pa- and p^-plaquettes belong to the same cube. It is clear that all the constraints 
are independent. Thus, the integral @ is equal (up to a possible constant factor) to the 
product of (^-functions (@) taken over all cubes c. This has to be substituted into ([!]): 



z=j udu P e^ u ^us( n u PA , n U PB 



(9) 



p c PaGc Pb&c 

Thus, we formulated the original model (|lj) in terms of the plaquette variables only. 

We notice now that this formulation has an interesting feature. Namely, @) is the partition 
function of a certain two-dimensional model. To see this, let us consider the section of a 
cube by hexagon, as shown at Figjl|. Using the one-to-one correspondence between the 
plaquettes of original lattice and the links of the hexagon, we replace the C/p-variables in (||) 
by (new) U[S where I is the corresponding link of the hexagon. 

The hexagons form the 2d lattice (Fig.|2|). The partition function (|9|), in terms of this 
lattice, takes the form 



z=f n^ tr (^n< n ^> n ^ 

/ h heh heh 



(10) 
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where I and h denote links and hexagons of the 2d lattice. The order of links in any hexagon 
(argument of (^-function) is hhhhhh (at Fig. 2(a), the h-type links correspond to the dotted 
lines). The fragment of the resulting 2d lattice is shown at Fig. 2(b). The whole lattice is 
easy to imagine - it is obtained from the regular 3d lattice (the lattice spacing is doubled 
with respect to the original 3d lattice) by replacing links by the tubes (handles) and vertices 
by the smooth connections of tubes. The surface separates all A-vertices from all £>-vertices. 
The two sets, being interior and exterior to the surface, are identical to each other. The 
genus g of the surface is proportional to A3, namely, g = gA/3. 

Thus, we started with the three-dimensional field and discovered that, actually, interaction 
takes place on a certain surface. In our view, the picture becomes more clear in terms of the 
r field (where r is the iV-component highest weight vector of representation; see definition 
(H)). Each particular sum over r is attached to the particular hexagon and we can consider 
(|l0|) as the partition function of statistical model of the iV-component 2d field r ((|i~0|) is the 
sum over all r-configurations) . 

It is clear from the above that the reduction described is a unique property of the gauge 
field and takes place due to the gauge invariance of the original model <^j. However, it is still 
interesting to see what happens in other models if similar change of variables is performed. 

We consider the principal chiral field (PCF) model. Its partition function has the form 

Zpcf= [HdUiUdUte^^+^SiU^i^UiU]) , (11) 

i I 

where i,j are the vertex labels, I denote links and I = {ij} denotes the link between two 
neighboring vertices i and j. 

Now our aim is to formulate this model in terms of link variables only, i.e., to perform 
integration over vertices. Due to the Euler theorem, A/"o — A/i + A/2 = 0, the topological 
integral over vertices produces A2 (^-functions (or A/2 — A3 + ... (^-functions, if dimension is 
higher than two). The corresponding condition for the plaquette is, obviously, 

s(i,U u i) ( 12 ) 

lep 

(this condition is valid for any loop, but we still keep the notation p). Thus, 

z PCF = I n w+B?) n ^(/, n ^) • w 

J 1 v ie P 

Here, p counts all independent loops, where 'independence' of the loop means an independence 
of the corresponding constraint. In particular, in the 2d PCF all the loops are independent 
(unless this 2d is the compact surface, in which case only one loop has to be excluded). 

Let us consider PCF on the surface shown at Fig. 2, in order to compare with QCD3. 
There are conditions 5^/, FJ/ U^j at the non-contractible loops (the non-trivial cycles). Field 
r corresponding to such a loop is attached to the section of the tube, i.e., there is a three- 
dimensional interaction inside the handles, in sharp contrast to QCD. 

If in (|l3|) we replace the <5-function by the heat kernel 

f(a,U) = J2d r e-^ r) X r(U) (14) 
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(C^r) is the second Casimir eigenvalue) and introduce a different couplings, a^ ex for the 
hexagons, and a tu b for the sections of the tubes, then Zpqf is restored as the Oih ex ,a. tu b — * 
limit of the result, while QCD3 partition function (|i~0| ) is the ah ex — » 0, a^j, — > 00 limit. 

There are several comments in order. The first one is the technical comment on use of 
(B) and (|5|) instead of (0). It has been understood already in QCD2 || that the O-degrees 
of freedom, though they are hidden in the model, play an important role. In QCD2, they 
act by means of the formula ^ which pulls the powers of d r to denominator and, thus, the 
1/iV-expansion of the QCD2 partition function || can be done ||. Here, in three dimensions, 
the use of the O-degrees of freedom allowed us to decrease the effective dimension of the 
system. For reader familiar with the hermitian matrix models we remind that the similar 
situation takes place there. In all these cases, the use of the ^-variables has led to the kind 
of a 'stringy' description. It is worth to remind that it is these variables, which are precisely 
absent in any abelian model. 

The important question is what happens with this picture in continuum limit. The loop 
variables are only non-zero observables in continuum theory as well. Therefore, the principal 
possibility to formulate the continuum model completely in terms of the loop variables still 
exists ||. From the present consideration we see that the change of variables from links to 
plaquettes can be done at any lattice scale. Suppose we approach the continuum limit by 
refinement of the 3d lattice so that the model is restored at each finer level. This should 
be achieved by the appropriate tuning of the coupling constant. Then at each level we can 
restore the described 2d picture with the coupling constant corresponding to the current 
level of refinement. However, the 2d lattice undergoes the infinite change of genus under 
refinement of the original 3d lattice. The quantity which survives the refinement is the ratio 
A/3/5 = 8/3, where A/3 is the volume of the original lattice. 

Note added: After submission of this work, I was informed of the related, pioneering work 
H on field strength and plaquette variable formulation of gauge theory. The relation between 
this early work and present paper is currently under investigation. 
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